Mesoscopic conductance fluctuations in graphene 
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We study fluctuations of the conductance of micron-sized graphene devices as a function of the 
Fermi energy and magnetic field. The fluctuations are studied in combination with analysis of weak 
localization which is determined by the same scattering mechanisms. It is shown that the variance 
of conductance fluctuations depends not only on inelastic scattering that controls dephasing but 
also on elastic scattering. In particular, contrary to its effect on weak localization, strong intervalley 
scattering suppresses conductance fluctuations in graphene. The correlation energy, however, is 
independent of the details of elastic scattering and can be used to determine the electron temperature 
of graphene structures. 

PACS numbers: 73.23.-b, 72.15.Rn, 73.43.Qt, 81.05.Uw 



INTRODUCTION 

Fluctuations of the conductance of electron systems 
have been studied for many years in metals and elec- 
tron gases formed in semiconductors They origi- 
nate from the interference between phase-coherent elec- 
tron paths that cross the system, and in diffusive, phase- 
coherent systems the conductance fluctuations are uni- 
versal (UCF): their average amplitude of the order of 
e 2 /h is independent of the mean conductance Qj. An- 
other consequence of quantum coherence in diffusive sys- 
tems is a weak localization (WL) correction to the con- 
ductance. It was shown recently that WL in a new two- 
dimensional system, graphene is unusual in that it 
depends not only on inelastic but also elastic scatter- 
ing of chiral carriers, both within and between the two 
graphene valleys 0, [|| . 

As quantum interference lies at the origin of both UCF 
and WL, it is important to understand if manifestation 
of UCF in graphene is also different from that in con- 
ventional two-dimensional systems. In our earlier stud- 
ies of WL in monolayer and bilayer graphene [H, 0] we 
have established that quantum interference of carriers 
within one valley is suppressed, but the presence of sig- 
nificant intervalley scattering makes it clearly detectable 
in the magnetoconductance. Here we perform the first 
analysis of conductance fluctuations in graphene, both in 
the metallic and electro-neutrality (Dirac) regions, which 
are complemented by WL studies. Using recent theories 
of UCF in graphene @, H, @|, we show that, similar to 
WL, the amplitude of the fluctuations can be strongly af- 
fected by the intensity of elastic scattering, so that it can 
be larger than in conventional 2D metals. In graphene 
flakes, however, significant intervalley scattering reduces 
it to a value close to that in one-valley 2D systems with 
non-chiral carriers. In order to detect the unusual prop- 
erties of UCF in graphene with a significant (up to four 
times) increase of the variance, one has to fabricate sam- 
ples with no intravalley suppression of interference as well 



as weak intervalley scattering. 

We also examine the effect of elastic scattering in 
graphene on the autocorrelation function of the fluctu- 
ations as a function of the Fermi energy, and show that 
under usual experimental conditions the correlation en- 
ergy is insensitive to the specifics of scattering and can 
be used as a direct measure of the electron tempera- 
ture. Methods to determine the true electron tempera- 
ture become important in graphene-based devices where 
phonons are poorly coupled to the environment, so that 
with commonly used currents the overheating of elec- 
trons is highly likely. The classical conductivity has a 
very weak temperature dependence and therefore is 
not suitable for this purpose. Both the WL and the vari- 
ance of UCF have strong temperature dependences [H, 0] , 
but the magnitudes of these effects depend on the details 
of inelastic and elastic scattering and so do not readily 
allow for the extraction of the temperature from them. 
We show here that analysing the width of the correlation 
function does unambiguously yield the electron temper- 
ature in graphene, as well as in any mesoscopic diffusive 
conductor. 



THEORY 

The variance of the conductance fluctuations in 
graphene is determined by the standard set of the per- 
turbation theory diagrams in h/(pp£) <§; 1 (£ is the 
electron mean free path) shown in Fig. [T] Structurally, 
they coincide with the diagrams describing mesoscopic 
fluctuations in usual conductors [lH but the Hikami 
boxes are different 0, 0, E3] because of the linear en- 
ergy spectrum, carrier chirality and valley degeneracy 
in graphene. In the low-temperature limit, Lt 3> L v , 
where Lt = {hD /hsT) 1 / 2 is the thermal diffusion length, 



L v = {Di 



a/2 



, (D is the diffusion coefficient and t v is 



the phase-breaking time), the variance of conductance 
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fluctuations becomes 0, HI 

(SG 2 ) = = ^5 s 2 ^(i, W, L v , L u L.) , (1) 

where G is the conductance in units of e 2 //i, 8G = G — 
(G), W and L are, respectively, the sample width and 
length, g s is the spin degeneracy, and the time-reversal 
symmetry parameter j3 — 1 for B -C Bq = (fi/e)/L 2 and 
/3 = 2 for fl > B Q . The function 7£(X, W, L p , Li, L*) is 
given by 
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( ?7 „„ i + 2(L/L l ) 2 )- 2 + (2) 
2(r] nm + ((L/Li) 2 + (L/L„) 2 ))- 2 } 

= + ^ 2 + 2^3 , 

where 77„ m = ir 2 n 2 + (L/W) 2 ir 2 m 2 + (L/L v ) 2 , Li = 
(Dt/) 1 ' 2 is the intervalley scattering length determined 
by the scattering rate t" 1 , and L* = (D / (t~ x + t~ 1 )) 1 / 2 
stands for the intravalley scattering length determined 
by both the trigonal- warping scattering rate r" 1 and the 
scattering rate r" 1 due to sublattice asymmetric poten- 
tials 0. 

In the case of L v 3> L the conductance variance given 
by Eqs. Q] and [2] is independent of the dephasing length 
and can be represented as 1Z = alZi, 4 > a > 1. The co- 
efficient a is sensitive to the strength of intervalley and 
intravalley scattering in a particular graphene sample. 
Namely, a = 4 if Li, L* 3> L; a = 1 for L ^ Li, and 
a « 2 for Li > L > L*. Therefore, the variance of 
conductance fluctuations in graphene is different from 
that in a conventional metal 2], in which case 1Z = 1Z\ 
in Eq. [5] and the variance is insensitive to the micro- 
scopic details of disorder. In the case of a long and nar- 
row graphene sample, strong intervalley scattering at the 
edges (L 3> Li) results in a = 1 and 1Z\ = Q(4)/i: A , 
leading to (SG 2 ) = (2//3)g 2 /15, which coincides with the 
usual result for quasi-lD disordered wires. 

At finite temperatures, the autocorrelation function 
of conductance fluctuations, F(A), is given by a con- 
volution of the ensemble-averaged correlator, J~(e) = 
(SG(E) SG(E + e)), and the thermal broadening factor 



K(e,A)= J &Ef(E,,. 



,)f'(E + e,e F + A), 



where f'(E,e F ) - -l/(4fc B T e ) cos\T 2 [(E - e F ) / (2k B T e )] 
is the energy derivative of the Fermi-Dirac distribution 
function (T e is the electron temperature). Thus 



F(A) = ((5G(e F )SG(e F + A)» 
= J deK(e, A)^(e) . 
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Here the brackets ((■■■)) stand for both the ensemble 
and thermal averaging and 
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T%* = [-jfTD + Vnrn + (L/Li) + (L/L*) 2 ) , (7) 

where td = L 2 /D. 

The variance of the conductance in graphene and 
the correlator T(e) given by Eqs. HE are sensitive to 
the microscopic details of disorder in a particular sam- 
ple. However, in the normalized autocorrelation function 
F n = F(A)/F(0) these details do not manifest them- 
selves under usual experimental conditions. 
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FIG. 1: (a),(b) The diagrams which contribute to the main 
order in the diagrammatic expansion of the conductivity- 
conductivity correlation function. The solid lines are the im- 
purity averaged retarded or advanced Green functions, the 
short wavy tails are the current vertices and the long wavy 
lines the diffusion ladders. (c),(d) Hikami boxes of two types 
and additional diagrams which determine renormalization in 
the main order in l/kpl <C 1, where the dashed lines corre- 
spond to the disorder potential, (e) Diagrammatic equation 
for renormalized current vertex. 



To show this, let us first consider a narrow wire, 
W -C min(L ¥ ,,L), with strong intervalley scattering at 
the boundaries. Then the part of the correlator T which 
contributes to F n can be written exclusively via the 
'valley-singlet' diffusion propagators T> 1 lm . The sum over 
77i in Eq. [4] in this case is dominated by the m = 
term. (Note that we have taken into account only the 
diffusion modes and neglected the Cooperons, which cor- 
responds to the regime of suppressed WL by magnetic 
field.) We obtain an asymptotic expression for F n (A) 
assuming L 3> L v keeping only the term with m = and 
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performing the summation over n: 
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where t = \J (eT v /h) 2 + 1. This function is sharply 
peaked, with a maximum at e = and width h/r v . In 
contrast, the thermal broadening factor K(e, A) has a 
broad peak at e — A of width ~ k^T e ^> S/r„. Convolu- 
tion of these functions results in a normalized correlation 
function that is independent of the microscopic details 
contained in !F(e): 



^n(A) = 



K(Q,A) _ 3(0coth0- 1) 
^(0,0) " 



A 



sinli e 



2k B T e 



(9) 



The width of this one-parameter function, determined 
from F n (9 c ) = 0.5, is 9 C = 1.36, which results in 



2.7fc B T e 



(10) 



This result remains valid under the condition Lx,L y <C 
min(L v , L x ) for any L v < L x , and allowing for all the 
diffusion modes in Eqs.[S][7]in graphene or, indeed, in any 
other mesoscopic disordered conductor. This result was 
tested numerically in 0. There, the values of A c were 
found to lie within a narrow interval, 2.7 < A c /k B T e < 
2.9, within 10% of the asymptotic value 2.7 of Eq. flOl 
Good agreement with this value (within 25%) was also 
found when 2D samples were considered. 



EXPERIMENT 

The samples studied experimentally in this work are 
monolayer graphene flakes created by mechanical exfo- 
liation on a n + Si substrate covered by 300 nm of SiO^- 
In addition to monolayer graphene flakes, we have also 
studied a bilayer sample (Bi: L = 1.5 /xm, W = 1.8 /im, 
see @), to examine the generality of the method of deter- 
mining the electron temperature. (The number of layers 
in the studied structures has been established from the 
analysis of the quantum Hall effect @, @j.) The manifes- 
tation of WL effects in monolayer and bilayer graphene 
can be very different, but the correlation properties of 
conductance fluctuations as a function of the Fermi en- 
ergy are expected to be universal. 

The ac current driven through the sample was 1 nA 
to avoid overheating (determined by measuring the ef- 
fect of increasing the current on the mesoscopic fluctua- 
tions). The Fermi energy, ef, was controlled by a gate 
voltage, V g , applied between the substrate and the flake 
(in our monolayers, Sf — frvp y/TrVgC/e, where C is the 
capacitance per unit area between the gate electrode and 
graphene). The conductance of all samples is shown as 
a function of Fermi energy in Fig. [5] The dephasing rate 
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FIG. 2: Dimensionless conductance of the graphene samples 
as a function of the Fermi energy at T =0.25 K. (For B2 the 
hole region is studied so the scale has been reversed to com- 
pare with other samples.) The shaded regions show the energy 
range used in the analysis. The inset shows the distribution 
function of 5G for sample F2 averaged over ~ 14000 'realisa- 
tions' (filled symbols are for high density, open symbols for 
the Dirac region). 
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0.4 0.06 
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F2 


3.8 


1.8 


3.8 


1.7 


0.5 0.1 


0.1 


0.7 



TABLE I: Lengthscales (in microns) of samples Fl and F2 
at high carrier density. The values of the dephasing length, 
L v have been determined from analysis of the weak localiza- 



tion. Shown here are values at T = 0.25 K (L" 



and 4K 



(ZA K ). i is the mean free path and n is the carrier density in 
10 cm 



was determined from a fit of the magnetic field depen- 
dence of the sample conductance to the theory of weak 
localization in graphene (Indeed, in order to anal- 

yse the UCF properly the values of not only L v but also 
Li and L* need to first be determined and this can only 
be done in combination with analysis of the WL.) The 
sample parameters are given in [9( and extra details of 
samples Fl and F2 are given in Table |TJ Samples Bl and 
B2 are quasi-lD samples with width W = 0.3 /jm and 
length L = 3.7 and 2.0 /zm, respectively. 

Fluctuations of the conductance occur as a function of 
both the Fermi energy and perpendicular magnetic field. 
The fingerprint of the fluctuations is robust and temper- 
ature dependent over the whole experimental range of 
0.26 to 20 K. For a particular range of energies and fields 
the autocorrelation function F(A) of the fingerprints can 
be analysed to determine the variance F(0) and the cor- 
relation energy A c from F(A C ) = F(0)/2. The V g range 
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was chosen to incorporate a sufficient number (~ 100) of 
fluctuations without significantly changing the average 
value of the conductance. In our experiments this limits 
the range of temperatures to T < 10 K. To increase the 
number of fluctuations (sample 'realisations') the mea- 
surements were performed at several (10 — 20) values of 
magnetic field in the range B 3> Bq. 




F n (A) 



FIG. 3: (a) Typical fingerprint of SG (normalized by e 2 /h) 
for sample Bl at T = 0.25 K and B — 90 mT as a function of 
Fermi energy, (b) Normalized correlation function for sample 
Bl. The inset shows the circuit used in the experiment and 
construction of the graphene sample. 



Figure Ufa) shows a typical fingerprint of conductance 
fluctuations (with average background, (G), removed) in 
sample Bl as a function of the Fermi energy in a magnetic 
field of 90 mT. In Fig. [3jb) the (normalized) autocorre- 
lation function of the fingerprint in (a) is shown. In the 
inset of the figure the circuit used in the measurements 
is shown. 

The random nature of the fluctuations is demonstrated 
for sample F2 in the inset of Fig. [2] where the distribution 
of the conductance fluctuations is seen to have a Gaussian 
shape typical of UCF [l3| ■ (It is also seen that for this 
sample the magnitude of the fluctuations is very similar 
at high carrier density and in the Dirac region.) Figured] 
shows the variance of samples Fl and F2 at high carrier 
density. One can see immediately that the variance is 
temperature dependent over the whole range of temper- 
atures down to 0.25 K. In the WL experiment, however, 
it was found that the magnitude of the quantum correc- 
tion determined by L v saturates at low temperatures [5| 
(when L v ~ L). This highlights the importance of the 
lengthscale Lt for the amplitude of UCF in 2D, when 
L T < fil. 




T(K) 

FIG. 4: The variance of the conductance fluctuations of sam- 
ple Fl (circles) with those of sample F2 (squares). Solid sym- 
bols correspond to the experimental results while open ones 
represent theoretical values calculated using Eqs.[3][7]for F(0). 



ANALYSIS AND DISCUSSION 



Using the values of L v , Li and L* determined from 
weak localization measurements in Eq. [3] we are able to 
calculate numerically the theoretical values of the vari- 



ance for our samples. For Fl and F2, L ¥ 



L > Li > £» 



One can see from Eq. [2] that the last two terms in TZ are 
suppressed by the strong elastic scattering and therefore 
a ~ 1. The good agreement between theory and ex- 
periment, Fig. 4, shows that L v calculated from WL 



and UCF is the same quantity, as predicted in [14J, and 
also that unlike its effect in WL intervalley scattering 
acts to suppress the conductance fluctuations. In || it 
was shown that there is strong intravalley suppression 
of quantum interference in graphene deposited on silica, 
hence the four-fold increase of conductance fluctuations 
(Theory section) will not be seen. In addition, due to 
intervalley scattering at the sample edges the interval- 
ley diffusion length Lj has an upper limit of the sample 
width. Therefore in narrow graphene samples it is diffi- 
cult to achieve conditions where a > 1, and as a result the 
amplitude of fluctuations is similar to that observed in a 
usual diffusive metal. We found previously [l6| that the 
fluctuations in bilayer sample Bi can also be interpreted 
using a standard UCF model. This shows a dramatic 
contrast with WL where in the presence of intervalley 
scattering the shape of the magnetoconductance curve 
shows strong dependence on small variations in the elas- 
tic scattering. 

The values of the correlation energy extracted from 
the correlation functions are shown in Fig. [5^a). It can 
be seen that the correlation energy increases linearly with 
increasing bath temperature. Equation [TU] is plotted in 
Fig. [5] as a solid line. There is agreement between the 
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FIG. 5: (a) Experimental value of the correlation function 
width A c for various bath temperatures extracted from the 
UCF in four graphene-based devices at high carrier density 
(F2: open squares, Bl: filled circles, Bi: open triangles). The 
solid line corresponds to Eq. 1101 (b) Correlation magnetic 
field of sample B2. Filled squares are experimental results 
and open squares are AB C — (h/e)/W L v . The dashed hor- 



izontal line shows AB C — (h/e)/W L. (c) Ratio L v /Lt 
high carrier densities in samples F2, Bl and Bi. 



for 



theory and experiment for all the samples. This includes 
the bilayer sample (open triangles in figure), which con- 
trasts with the significant difference between monolayer 
and bilayer samples in the manifestation of weak local- 



ization 



ndicates general applicabil- 
ity of the developed model to all diffusive systems. The 
criterion of the applicability of the theoretical result in 
Eq. [lOl L^/Lt > 1, is seen to be satisfied in the studied 
samples, Fig. EJc). 

We want to note that the results of our experiments on 
graphene show some general properties of UCF |2|. In- 
creasing magnetic field above Bq has shown the conven- 
tional decrease of the variance by a factor of two. Also, 
the correlation magnetic field has been seen to be related 
to L V (T) in the usual way. Shown in Fig.^b) is the cor- 
relation field for sample B2 as a function of temperature. 
For this sample L v > W, so the expected dependence 
of the correlation field is AB C (T) » (h/e)/WL v (T) fl. 
This is plotted as open squares in the figure and clear 
agreement is seen with experiment, including the satura- 
tion of AB C that occurs at low temperatures when the 
dephasing length approaches the sample length. 

Finally we discuss the results of our measurements in 
the Dirac region where the applicability of the developed 
analytical theory of UCF is not obvious, although the 
conductivity of our samples satisfies the commonly used 
criterion for the diffusion theory, (G)L/W 3> 1. (We 
have also shown previously || that the diffusive theory 
of weak localization describes well the magnetoconduc- 
tance in the Dirac region.) Numerical investigation of 
the conductance fluctuations in graphene (from sample 



to sample) [19( have shown that their amplitude is consid- 
erably stronger than in conventional metals. At the same 
time experiments on few-layer samples [2(| has claimed 
that conductance fluctuations (as a function of the car- 
rier density) are suppressed near the charge neutrality 
point. 

In our experiments we consider a region of carrier den- 
sities around the Dirac point, —30 < ep < 30meV. 
Firstly, we do not observe an increase of the amplitude of 
conductance fluctuations in the Dirac region compared 
with that in the high-density regions. In samples Fl 
and F2 the variance in the Dirac region is very close 
to that in high-density regions, while in narrow samples 
Bl and B2 it is several times smaller. Secondly, it can 
be seen from Fig. [6] that the correlation energy deviates 
from the 'high temperature' limit at low temperatures. 
This can be attributed to the fact that for all samples 
the high-temperature condition is now destroyed because 
of a decrease of the dephasing length by > 30% that 
occurs when moving from high to low carrier density, 
Fig. [Bpnset). (This decrease of the low-temperature de- 
phasing length in the Dirac region was obtained from the 
analysis of WL in both monolayer and bilayer graphene 
samples @,E3]-) 
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FIG. 6: A c as a function of bath temperature at carrier den- 
sity close to the Dirac point (F2: open squares, Bl: filled 
circles, Bi: open triangles). The solid line corresponds to 
Eg. 1101 Inset: Ratio L v /Lt at low carrier density in samples 
F2, Bl and Bi. 



CONCLUSION 

In conclusion, in combined studies of conductance fluc- 
tuations and weak localization in monolayer and bilayer 
flakes we have demonstrated that the variance of univer- 
sal conductance fluctuations in graphene is strongly af- 
fected by elastic scattering, in particular intervalley scat- 
tering. However, the correlation energy of the fluctu- 
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ations as a function of the Fermi energy is insensitive 
to these scattering mechanisms under common experi- 
mental conditions. Analysis of this correlation energy al- 
lows direct measurement of the electron temperature in 
graphene structures. We have also discussed the evolu- 
tion of the variance and the correlation energy of the con- 
ductance fluctuations from the region of high carrier den- 
sity to the region close to the point of electro-neutrality. 
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